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WE CLAIM: 

1. A method for the generation of three dimensional fractal subsurface structure by Voronoi 
tessellation and computation of gravity response of such fractal structure, by generation 
of fractal subsurface structures which are very close to the natural setting of the 
subsurface geology and provides computation of forward gravity response of such 
structures for delineation of the underlying anomalous object, said method comprising: 

(a) selecting Voronoi centers at a plurality of locations over a region of interest, the 
Voronoi centers being represented by X, Y, Z co-ordinates; 

(b) generating an initial model of the subsurface fractal geological object, having 
variation in physical property such as density in lateral and vertical directions; the 
initial model being generated by tessellating, the Voronoi centers, and assigning 
values of physical property variations during generation of the model on the basis of 
pre-determined assumptions in the model; 

(c) assigning different colours to regions in the model which have different physical 
properties such as density to enable demarcation. 

2. A method as claimed in claim 1 wherein the fractal subsurface model is generated by 
modified Voronoi tessellation technique which comprises modifying the Voronoi 
tessellation by using L p norm, where p is an exponent which can assume any real value, 
thereby enabling greater range for generation of different fractal sub-surfaces, which are 
very close to the real geological situations. 

3. A method as claimed in claim 1 wherein the natural setting of the geological subsurface is 
selected from the group consisting of a sedimentary basin, hydrocarbon deposits, oil 
reservoirs, aquifers and mineral deposits. 

4. A method as claimed in claim 1 wherein the computation of gravity response due to 
fractal subsurface generated by modified Voronoi tesselation is carried out by analytical 
expression which comprises: 

(a) demarcating boundaries of regions having different physical property, such as density 
in the tessellated region, the boundaries forming a polygonal shape in 2-dimensional 
space; 
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Where S= +1 if pi is positive, S= -1 if pi is negative, 

W= +1 if mi is positive, W= -1 if mj is negative, 

'Z' is depth and V is number of sides in the polygon. 

G is universal gravitational constant, p is the density of the tessellated regions; 
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where the effect of common arm of the adjacent polygon is removed; 

(b) repeating the gravity response computation due to the polygon for all adjacent 
polygons of different physical properties using the demarcated polygon boundary; 

(c) adopting the process of gravity response computation for tessellated regions lying at 
different depths to obtain a response; 

(d) integrating the response using Simpson/Gauss quadrature formula at plurality of grid 
nodes overlain on the region of interest. 

A method as claimed in claim 1 wherein the results are achieved through self-written 
software, which generates Voronoi tessellated subsurface region and computes gravity 
response of the same. 

A method as claimed in claim 1 being applied to geophysical inversion, wherein the 
tessellated regions are altered by changing the position of Voronoi centers in each 
iteration. 
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